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We argue that WIMP dark matter can annihilate via long-lived "WIMPonium" bound states in 
reasonable particle physics models of dark matter (DM). WIMPonium bound states can occur at or 
near threshold leading to substantial enhancements in the DM annihilation rate, closely related to 
the Sommerfeld effect. Large "boost factor" amplifications in the annihilation rate can thus occur 
without large density enhancements, possibly preferring colder less dense objects such as dwarf 
galaxies as locations for indirect DM searches. The radiative capture to and transitions among 
the WIMPonium states generically lead to a rich energy spectrum of annihilation products, with 
many distinct lines possible in the case of 2-body decays to 77 or 7Z final states. The existence of 
multiple radiative capture modes further enhances the total annihilation rate, and the detection of 
the lines would give direct over-determined information on the nature and self-interactions of the 
DM particles. 
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I. INTRODUCTION 



One of the most promising ways in which to probe the 
nature of weakly-interacting-massive-particle (WIMP) 
dark matter is provided by indirect detection experi- 
ments in which the annihilation products of dark matter 
in astrophysical contexts are observed. Indirect detec- 
tion experiments such as AMS [1], ATIC [2], EGRET 
[3], GLAST [4], HEAT [5], HESS [6], INTEGRAL [7], 
PAMELA [8], and VERITAS [9] look for signals ranging 
across energetic gamma rays, positron excesses, and anti- 
proton fluxes, and hints of deviations from background 
expectations now abound. A common feature in the in- 
terpretation of these experiments is that the WIMP an- 
nihilation proceeds via simple, almost free-particle anni- 
hilation leading to a rate that depends on the WIMP rel- 
ative velocity in only a very simple, essentially structure- 
less way, and moreover, is directly related to the cross- 
section that led to the DM density at freeze-out. This 
assumption then implies that the primary astrophysical 
quantity determining the fluxes from DM annihilation is 
the local DM squared-density, p(x) 2 , with the velocity 
distribution of the DM being essentially irrelevant. In 
addition the assumption of almost free-particle annihila- 
tion gives a relatively simple resulting energy spectrum 
for the annihilation products (though of course features 
such as steep falls at kinematic thresholds are generic). 
The dependence of the annihilation fluxes on p(x) 2 has 
been widely employed in the suggested interpretations of 
the various experimental anomalies as, very frequently, 



a so-called "boost-factor" in the annihilation rate is re- 
quired to match the observed flux, and this is assumed 
to come from local over-densities of DM. 

In this paper we point out that annihilation of WIMP 
dark matter via intermediate long-lived "WIMPonium" 
bound states, Q( n ' e \ is possible in many particle physics 
models of DM [10] (see [11] for another recent discussion 
of WIMPonium). As we argue below, the WIMPonium 
bound states can occur at or near threshold in which case 
they lead to potentially very substantial (factors of 10 3 to 
10 5 ) enhancements in the DM annihilation rate, closely 
related to the well-known Sommerfeld non-perturbativc 
enhancement [12] that has been applied to freeze-out 
calculations and indirect dark matter signals [11, 13- 
16]. The existence of this p 2 - independent dynamically- 
induced "boost factor" has important implications: First 
large amplifications can occur in low-velocity-dispersion 
systems without a large (eg, cusp-like) density enhance- 
ment , possibly preferring colder less dense objects such 
as dwarf galaxies as more promising places to search for 
DM signals compared to the higher density but higher 
velocity galactic center 2 . Second, a knowledge of the 
full phase space density of the DM is necessary to reli- 
ably compute the DM annihilation rate, which strength- 
ens further the case for detailed realistic simulations of 
the DM distribution in our galaxy. 

In addition, the deeply-bound WIMPonium spectrum 
is often quite involved, and the radiative capture to and 
transitions among the various states generically lead to 
a rich energy spectrum of annihilation products, with 
many distinct lines (in the case of decays to 77 or "fZ 
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1 Furthermore cusp-like dark matter distributions seem not to be 
favoured at present, see for example [17]. 

2 See [18] for a travel guide to the directions in the sky where DM 
annihilation signals may be expected 
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final states, similar to that found in WIMP annihilations 
[19]) possible. Although the observability, or otherwise of 
these distinct lines depends on the DM model being stud- 
ied, and especially the resolution of the detector, the ex- 
istence of all the various radiative capture modes further 
enhances the total annihilation rate, and the detection of 
even just a few of the lines would give direct information 
on the interactions of the DM particles. 



II. THRESHOLD BOUND STATE BASICS 

In order to demonstrate the important features of 
WIMPonium bound states it is useful to examine a sim- 
ple, almost model independent, example. We introduce 
two complex scalar fields s and n with the following in- 
teractions and masses 



X'fj, .ml 2 ml 
C = —nss + — \n\ + — 



(1) 



where A' is a dimensionless coupling and /x, m n and m s 
are mass dimension one parameters. We assume that 
e = m n /m s <C 1 (typically we will take m n ~ mz, while 
500 GeV < m s < 30 TeV) . A discrete symmetry is im- 
posed forcing s particles to appear in the Lagrangian in 
pairs so the s scalar will be the stable dark matter WIMP. 

As we have indicated in the introduction we are inter- 
ested in two related phenomena: 1) The scattering of slow 
moving particles near bound state thresholds leading to 
amplification of the direct annihilation rate; 2) The ra- 
diative capture to and decay of deeply bound states and 
transitions between different bound states. 

Since we are interested in low-velocity processes we 
can proceed by solving the Schrodinger equation for 
the two-s-particle system with the Yukawa potential 
V = —\ 2 e~ mnr /8irr that follows from the exchange of 
n-scalars. Here A = X'fi/m s . Semi-classical considera- 
tions show that the number of bound states Ne of given 
orbital angular momentum I satisfies [20] (21 + l)N e < 
2M r J r\V(r)\dr, where M r is the reduced mass of the 
two particle system, so, for example, the number of S- 
wave bound states in our case satisfies N < ^-a, where 
a = \ 2 /8ir. A more precise condition for there to be 
at least one bound state follows from numerical methods 
giving [21] a > 0.84m„/m s . Thus if we want to have a 
rich structure of energy levels our dark matter particles 
will cither need to have large couplings with the n scalar 
or need to have m s /m n >• 1, or both. We emphasize, 
however, that the most important phenomenology - the 
large amplification of the DM annihilation rate - requires 
only a single at-threshold bound state, and so imposes 
only a mild condition on the coupling. For instance for 
m n — mz, m s = 500 GeV, we require a > 0.15 which is 
well within the perturbative regime a < 2tt. Note that, 
upon writing the complex scalar s in terms of its CP- 
even and odd parts s — 4> s + ia s , we have, due to Bose 
symmetry, that the 4> s 4> s and a s a s bound states can only 
possess even orbital angular momentum I = 0, 2, while 



the angular momentum of the (f> s a s bound states is un- 
restricted. 

Turning to the scattering of two slow moving s par- 
ticles by the Yukawa potential arising from the n scalar 
exchange interaction, both elastic and inelastic scatter- 
ing cross sections (such as radiative capture) are ampli- 
fied by a non-perturbative Sommcrfcld-like enhancement. 
This enhancement can be formulated in terms of a non- 
relativistic quantum two-body problem with a potential 
acting between the incoming particles. To a good approx- 
imation this leads to a dressing of the dominant S'-wave 
part of the tree-level cross sections by a multiplicative 
factor, 



n 2=0 
^ u tree ' 



(2) 



An exact analytic calculation of R for a Yukawa poten- 
tial is not possible (although we give a close analytic 
approximation later) and so we must proceed numeri- 
cally. The Schrodinger equation for the radial part of the 
two dark matter particle state, ip(r), with I = 0, reads 
—ip"(r)/m s + V(r)ip(r) = Eip(r), where E = m s (3 2 is 
the kinetic energy of the two dark matter particles in the 
center-of-mass frame, where each dark matter particle 
has velocity 0. Using the outgoing boundary conditions, 
?A'(oo)/V>(oo) = im s /3 7 R is given as R — \ip(0) /ip(oo)\ 2 . 
Rewriting r — y/m n and letting e = m n /m s we can 
rewrite the Schrodinger equation as a function of the two 
ratios a/e and a//3, viz 



dr 2 



a e 



V 



(3) 



The resulting numerical solutions for R are functions of 
t = a/e and u = a/ f3 with the 2d contour plot shown in 
Fig 1 and the 3d version in Fig 2. 

It is clear that there are two distinct regions in Fig 1. 
For large values of the velocity (/3 > e) there is a relatively 
flat region - the Coulomb region. This is the part of 
parameter space relevant for freeze-out, the magnitude 
of the enhancements being at most a factor of 3 to 5 
[13, 14]. More interesting is the low- velocity region in 
which we see the effect of resonance peaks. These peaks 
correspond to the formation of I = bound states at 
threshold (E = 0), and are the focus of this analysis. 

Provided we are sufficiently close to a resonance peak, 
the dependence of R on t = a/e and u = a/f3 is described 
by a modified Breit-Wigner resonance formula applicable 
for threshold resonances due to Bethe and Placzek (BP) 
[22, 23]. Considering only elastic scattering, the Breit- 
Wigner resonance cross section is 



BW 



k 2 ((m s /3 2 -£ ) 2 +r 2 /4)' 



(4) 



where T e and T are the elastic dissociation and total 
width for the resonant bound state and So is the dis- 
tance of the resonance from exact zero-energy and is in- 
dependent of j3. Following BP, for small j3 near a thresh- 
old resonance the BW expression is modified by replac- 
ing T e = VE^e with 7 e not depending on E . Using 
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Contours in R 




LogWjS] derived from the numerical solution of Eq.(3). 



FIG. 1: Contour plot of the enhancement factor R. Shown is 
a typical "path" in parameter space as the velocity of the s 
states decreases from the value at freeze-out to that relevant 
for indirect detection. The path is not exactly horizontal due 
to thermal corrections to the masses and couplings. 
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FIG. 2: A 3d version of Figure 1. 



E = k 2 /m s , we have 



BP 



m s ((m s (3 2 - £ ) 2 + m s /3 2 j 2 /A) 



(5) 



(as R has been found by including only elastic scattering 
we have set T = T e ). 

Eq.(5) shows that for eo ^ the cross-section first 
increases as 1//3 2 and then becomes independent of (3 
for (3 C 1 as shown in Fig 3. The plateau begins at 
P ~ \/ko|/ m s, and corresponds to af p ~ 47r/(m s |£o|). 

Comparing the numerical calculation of the elastic 
cross section to a^ p in the low-velocity limit we can ex- 
tract the numerical values of j e and £q as a function of 



Resonance 


a 


to 


Is 


12.2 


1.681 


2s 


198.5 


6.453 


3s 


1015.9 


14.358 


4s 


3220.6 


25.407 


5s 


7849.0 


39.609 



TABLE I: Numerical fits of R(t) = a/(t - t ) 2 . The cross 
sections are extremely sensitive to the values of to- The values 
given have been rounded off so that they can be displayed. 



the parameters a and e for each of the possible near- 
threshold bound state resonances. 

The plateau arises as (3 — > (for Sq ^ 0) with asymp- 
totic value of the cross section 



1/3- 



7T 

TO., 



(6) 



Since a 



R(t)<r l tr ® e , where R(t) is independent of (3, 



7e 



-R(t)a tr 



(7) 



The function R(t) has the form R(t) = a^ nl) /{t - t K ™') 



= a (ni)/a + (n^)^2 

where and depend upon the principal and or- 
bital angular momentum quantum numbers (nl) of the 
resonance that is close to threshold as t is varied. Table 
I gives numerical fits for the S'-wave resonances. 

We remark in passing that an exact analytic treat- 
ment is possible for the Hulthen potential Vn(r) — 
— cvm„e _rm "/(l — e~ rm ") which has similar r — > and 
r — > oo behaviour to the Yukawa potential. The S-wave 
phase shifts are (see also [24]) 



7T 

2 



are 



r i 



ik 



l A 



(8) 
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where k = m s /3 is the momentum of the scattered state 
and A can be thought of as t — a/e up to a factor of two. 
From this the S-wave elastic cross section a e — 47r|/| 2 
follows using / = (e 2t5 ° - l) /2ik. Taking the (3 -> 
limit of a e , the dominant behaviour is a cx 1/(A — Ao) 2 , 
where A plays the same role as tg™ ^. Exactly on one 
of the resonances, A = Ao, the dependence on f3 is a e cx 
1//3 2 which agrees with the BP form. 

So far we have only discussed the case of elastic scatter- 
ing. However, as we are interested in the indirect signals 
coming from dark matter annihilations we need to ex- 
amine the case of inelastic scattering. In order to write 
down the inelastic cross section in the presence of the 
long range enhancements we again follow BP by writing 



BP 



T e T,, 



k 2 {{E~e a ) 2 + Y 2 /AY 



(9) 



where Ti is the inelastic width associated with the direct 
annihilation or radiative capture of the incoming s-pair. 
This form is again only applicable when we are suitably 
close to one of the S-wave resonances. Following BP the 
inelastic width is a constant as opposed to the energy 
dependent elastic width. Substituting the form for the 
elastic width T e = \J~E^ e into Eq.(9) we have 



BP 



pmV 2 [{m s (i 2 - e ) 2 + 



(r,+/3ym7 7e ) 2 



(10) 



The first point to note here is the 1//3 enhancement 
of the inelastic cross section compared to the elastic - 
the usual Bcthe 1/v law [25] . This follows if we take the 
j3 — > limit of Eq.(10) (assuming that we are not exactly 
on resonance) 



BP 



7. 



■Ti 



/3m 



l /2 (4 



r?/4)' 



(ii) 



compared to the f3 — > elastic cross section, erf p = 
7T7g/m s (£g +r?/4). Second, for m s j3 2 > s a ,Ti, the in- 
elastic cross section Eq.(10) shows a 1//3 3 dependence 
which plateaus to the standard 1//3 dependence when 
m s j3 2 <C So, Ti. This behaviour is exactly that of Eq.(2) 
for the enhancement of the naive inelastic cross section 
given the /3 dependence of R plotted in Fig. 3 resulting 
from numerical solutions and as discussed above. Third, 
the final limiting value of j3af p has a 1/e 2 , dependence 
and shows that for WIMPonium bound states close to 
threshold the cross sections are further enhanced. (In the 
above we have made the assumption that e > r» which 
is true over the majority of parameter space. However, 
for the case where £o < Ti the true dependences can be 
more complicated involving cross terms of the inelastic 
and elastic widths.) Depending on the size of e the size 
of the 1//3 2 enhancements can be significant since for 
DM annihilations in present day astrophysical systems 
the relevant range of (3 is <~ 1CP 3 — 10~ 5 . 

A useful way of thinking about and calculating approx- 
imately the amplification of the elastic and inelastic cross 



sections is in terms of the diverging scattering lengths 
that occur when a bound state energy tends to zero. Re- 
call from the elementary theory of non-relativistic elastic 
scattering that the S-wave phase shift satisfies 



lim fccot[f> (fc)] 

k— i-0 



1 1 

X + 2 



r e k 2 + ... 



(12) 



where r e <~ l/m„ is the effective range of the potential 
and l s is the scattering length which is related to the 
near-threshold bound state energy e by 



4 = 



l 



+ .... 



(13) 



These equations imply cot So = — tJ\so\/E where E — 
k 2 /m s is the CM scattering energy, and we have assumed 
k — > before £o — > 0. From the standard expression 
of the elastic cross-section in terms of So one then finds 
agreement with the BP formula Eq.(5) in the same limit 
of k — > before £o — > 0. Comparing with the BP form we 
learn that j 2 = 4|e | — > as the bound state approaches 
exactly the zero-energy threshold. The advantage of this 
approach is that the S-wave component of the distorted 
incoming plane waves can be simply expressed in terms 
of So and thus e - To a good approximation the wave- 
function is 



1=0 



sm(kr + So) 
kr 



(14) 



This form and its dependence on the scattering length £ s 
will be useful to us when we discuss radiative capture. 



III. DECAYS, CAPTURES, & TRANSITIONS 

WIMPonium bound states possess a rich phenomenol- 
ogy of radiative captures to various bound states, and 
decays and transitions from or among the bound states. 

A. Radiative Capture 

For interesting bound-state to bound-state transitions 
to be relevant, the radiative capture cross section must 
be significant. Similar to elastic scattering the radiative 
capture cross section is enhanced when there is a near- 
threshold bound state, as is approximately the situation 
in neutron-proton scattering where enhanced radiative 
capture to a bound deuteron is possible. 

We will consider transitions into both the £ = 
and £ = 1 bound state energy levels. The most eco- 
nomical way to perform such a calculation is to first 
expand the continuum state in terms of partial waves 
where we will keep only the S-partial wave, Eq.(14) 
(the incoming P-wave gives terms that are suppressed 
by /3 2 ). Labeling the bound state wavefunctions into 
which we will be capturing as ip n u z , the two most im- 
portant states are Vioo ~ ( 7ra o) _1 ex P(~ r / a o) and 
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V'zio ~ (327raQ) _1//2 r cos6'exp(— r/2ao)/ao, where we've 
assumed that the Vt^ n ^ bound state wavefunctions are 
similar to hydrogen with a ~ 2/m s a. (The bound 
states ^200 and ip2i±i do contribute, with the total rates 
for radiative capture changing by, at most, an 0(1) fac- 
tor. As we are interested in the general parametric de- 
pendence we neglect the contributions from capture into 
these bound states.) The radiative capture cross section 
depends on matrix elements of the form / = (ipf\0\ipi), 
where tpi(r) is the initial distorted partial wave of the 
continuum state, tpf(r) is the wavefunction of the bound 
state into which we are being captured, and is the in- 
teraction Hamiltonian. In our case — |e jp r , where p 
is the momentum of the radiated scalar n state. 

Expanding the exponential in powers of p.r and con- 
sidering radiative capture into the Is and 2p states, the 
overlap integrals take the forms 



J 2 - 

1 ls — 



A 2 

4 
\ 2 

j2 _ A 



J [V>ioo(pi-r) 2 V*=o] d 3 r 
J [^2io(P2-r)^=o] d 3 r 



(15) 



where the first non-zero integral comes at second order 
in p.r for capture into the Is state and at linear order 
for capture into the 2p state. 

From this the rates are found to be 



r 



2p 



288tt am 3 s (a oPl f(4a Q + S /k) 2 , 
1024tt am 3 s (a p 2 ) 3 (8a + 6 /k) 2 , 



(16) 



where p\ « \/E\ — and p 2 « \/E 2 — m 2 are the mo- 
menta of the n scalars due to transitions into the Is and 
2p states respectively. For simplicity, from now on we 
will assume that m n is small relative to transition ener- 
gies and so can be neglected. Clearly if m n is not small 
then there are trivial kincmatical suppression factors. 

Using the S-wave phase shift, Eq.(12), and taking the 
small k limit gives So ~ k£ s . For large scattering length 
(compared with 8ao) the rates become 



^2p 



97ra 6 m 3 J 2 w 9na 6 



2na 4 m 3 J 2 



27ra 4 



£o 
L rr^ 

£o 



(17) 



where in the last expression we have used the relation, 
Eq.(13), between the scattering length £ s and the near 
threshold bound state energy and continuing the analogy 
with hydrogen we have taken the bound state energies to 
be E n = —m s a 2 /4n 2 . This implies radiative capture 
cross sections 



RC 



RC 
°2p 



9na b 
/3s m s ' 

2na 4 
f3e m s ' 



(18) 
(19) 



This shows the usual 1//3 Bethe-dependence of an inelas- 
tic cross section near /3 — > 0, and most importantly an 
additional a 2 m s /so enhancement relative to the radiative 
capture cross section if there was not a near threshold 
bound state (in other words if the do- dependent terms 
dominated over So/k in Eq.(16)). The factor a 2 m s /e 
is just the ratio of the typical bound state energy com- 
pared to the near-threshold energy. In addition to this, 
it should be noted that for a ~ 0(1) (which is still well 
within the perturbative regime) the radiative capture 
cross sections are further increased, although the sim- 
ple hydrogen-like scaling that we have employed starts 
to break down and numerical methods must be used. 



B. Decays 

Consider annihilation of the ss bound state Q^ n '^ to 
light (mass -c m s ) degrees of freedom. Let the ampli- 
tude for the free 2 — > 2 scattering be Afi, then from the 
standard theory of, eg., positronium decay the amplitude 
Mfi for the bound state decay is 



M f 



1 

4m s 



3 



7T' 



/ 



d 3 p^(p)A fl ( P ), (20) 



where the momentum-space bound state Schrodingcr 
wavefunction is normalized as J d 3 p\ip{p)\ 2 — 1. Here pis 
defined by (0, 2p) = q 1 — q 2 where q± are the 4- momenta 
of the two s constituents. In the limit of relatively weak 
binding (Eb « m s ) we have Mq w 2m s . Expanding 
Afi in powers of p 2 , Afi(p) — A^ +p 2 A^ + ... and using 
i>(p) = J d z xe lp - x il){x)/{2T:f/ 2 gives 



M fi 



4°^(o) 



15^ & m 



dr 2 



(21) 



which enables the calculation of the decay width of vari- 
ous orbital angular momentum bound states. 

For m s /m n 3> 1, and the ssn coupling A not very 
large, the wavefunctions of the bound states are approxi- 
mately hydrogen-like, which upon applying Eq.(21) gives 
the annihilation width £ = 



£=0 



3n 3 



(22) 



This expression is most accurate for the more tightly 
bound WIMPonium states with small principal quantum 
number, n = 1,2,..., while for higher bound states a 
precise decay width requires numerical evaluation of the 
bound state wavefunction. Parametrically Eq.(22) gives 
a good estimate of the decay width in all cases. 

From Eq.(21) and the expansion of the free 2 — > 2 scat- 
tering amplitude, the decay widths of the higher orbital 
angular momentum states are parametrically suppressed 
by powers of 1/ (m s ao) 2 ■ Using hydrogen-like scaling once 



=o 

ann 



a which, if we take 



again shows that, eg, T^/T l a 
a < 1, is parametrically small. However, if we take 
a ~ 0(1) all decay widths can be large and compara- 
ble. 
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C. Transitions 

Transitions between the various bound states are pos- 
sible with either the emission of the CP = ±1 compo- 
nents of n, or if on-shell n-production is kinematically 
disallowed, by decay to light SM states through virtual 
n-emission. First assuming on-shell n emission (possible 
when a 2 m s > m n ), the relevant WIMPonium matrix el- 
ement reduces to T fi = (Rn>e>Ye> m >\ Aexp(ipr) \R n eY im ). 
Similar to transitions in hydrogen-like systems we may 
expand the exponential in powers kr ~ a (for relativistic 
n). The first transitions occur at 0(kr) with the emission 
of the CP — — 1 state of n with orbital angular momen- 
tum 1 = 1, changing the bound state from <j> s (f> s or a s a s 
to 4> s a s or vice versa. 

A good estimate of the various transition rates follows 
from a straightforward application of Fermi's Golden 
Rule. Parametrically the rate for At = 1 transitions 
emitting a relativistic a n scales as 

r^s 1 - aal (AE) 3 ~ a 5 m s (23) 

where AE (^> m n ) is the energy splitting between the 
bound states. We see the (AE) 3 behaviour familiar from 
hydrogen-like systems which favours deep transitions. All 
other transitions emitting an on-shell n are down by pow- 
ers of (AEao) 2 ~ a 2 . For instance Al = or 2 transi- 
tions emitting a CP = +1 <p n state are suppressed as 

r^fns ' 2 /^trans 1 ~ °^ as we ^ as ^ ^ na ^ state phase space 
factors (AE/AE') 3 . Transitions to light SM states via 
virtual n's are further suppressed by both couplings such 
as a or a em and (at least) three-body phase space factors. 

The most important feature of Eq. (23) is that it shows 
that the transition rate between different bound states 
can be competitive to the direct decay rate, Eq.(22), of 
the I = tower as long as on-shell n production is kine- 
matically possible. If the s DM particles are captured 
in a P-wave orbital angular momentum state, then the 
extra suppression of P-wavc annihilations implies that 
transitions via on-shall n's can dominate. Furthermore 
transitions to light SM states via virtual n's might be 
only mildly suppressed relative to P-wave annihilations 
depending on the strength of the coupling of n to Higgs 
(and thus other SM) fields. Fig. 4 schematically illustrates 
the dominant decays and transitions and their respective 
rates while Fig. 5 depicts diagrams responsible for discrete 
7 lines. 

IV. CONSEQUENCES FOR INDIRECT 
DETECTION AND LHC 

A wide range of consequences follow from the exis- 
tence of WIMPonium bound states. As we have already 
mentioned they provide a new, dynamical mechanism for 
the "boost factors" that are often introduced to explain 
anomalies in indirect detection observations. Unlike tra- 
ditional p(x) 2 enhancements they depend on the veloc- 
ity distribution of the DM particles, and since between 



1 = 1 = 1 1 = 2 




FIG. 4: Diagram illustrating, for various bound states, the 
rates of annihilation (double line arrows), and transitions 
emitting a 4> n or a n (single line arrows). We have omitted 
the 4> s a s I — and I — 2 bound state energy levels and the 
associated transitions and decays for clarity. 



oW) 




FIG. 5: Schematic diagrams showing annihilation of, or tran- 
sition between, different bound states leading to discrete 77 
or 7Z lines. The states in the triangle can be any charged 
scalar, W gauge boson, or fermion that couples to (f> n ,a n , 
in general via mixing with Higgs states. 



the freeze-out epoch and today the velocity changes from 
f3f <~ 1/5 to (3 now <~ 10~ 3 — 10~ 5 , the near-threshold 
bound state enhancement decouples the value of the an- 
nihilation cross section determined by successful thermal 
freeze-out from that observed today in indirect annihila- 
tion observations. (See [14] for a discussion of thermal 
freeze-out production of DM in models closely related to 
our toy theory.) In fact, because of thermal corrections 
to the couplings and masses during freeze-out the "path" 
taken in (a//3, oe/e) parameter space as the universe cools 
is not exactly a a/e ^constant line, but instead can move 
towards or away from the value at which a resonance 
occurs exactly at threshold. This further decouples the 
value of the cross-section at freeze-out from that observed 
now. Moreover, the enhancement of the cross-section as 
the universe cools leads to the possibility that there is in- 
teresting residual post-freeze-out annihilation of the DM 
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particles, for instance leading to changes in BBN predic- 
tions of such elements as 6 Li and 7 Li [26]. Depending on 
the maximum size of the threshold enhancement, which 
is determined both by the size of the inelastic widths, 
and the degree to which the bound state approaches zero 
energy, this dynamical boost factor can be large enough 
to potentially favour environments such as dwarf galax- 
ies, since their velocity dispersions go down as low as 
3-5kms _1 [27] compared to the typical galactic value of 
~ 200 - SOOkms^ 1 . 

Turning to particle physics model-building issues, the 
existence of WIMPonium bound states implies that the 
standard supersymmctric ncutralino DM picture must 
be modified somewhat. Although we have explained the 
phenomenology of WIMPonium in the context of a very 
simple, purely scalar toy model, we emphasize that simi- 
lar phenomena are possible if the DM is fermionic, or even 
neutralino DM. From the condition for the existence of 
at least one bound state a > 0.84to„/to s we learn that 
if the DM is a neutralino interacting via the exchange 
of and Z gauge bosons (and Higgs states) then the 
neutralino must be heavy m n cutraiino ^ mw jo-i ~ 2TeV. 
On the other hand if the DM particle interacts with a new 
strong-interacting sector, say a hidden valley sector[28], 
then the DM particle can potentially have close-to- weak- 
scale mass. 

A particularly attractive possibility is to have the DM 
state associated with electro-weak symmetry breaking 
dynamics in some way, so that it interacts with the Stan- 
dard Model via the so-called Higgs Portal [14], and has 
strong Higgs-mediated self interactions. This is in the 
same class of models as our toy theory, although the DM 
particle s could be fermionic in which case the spectrum 
of bound states and associated decays and transitions is 
even richer. In all three cases the LHC search strategy 
for the DM state is greatly modified compared to the 
standard expectation of weak-scale neutralino DM. 



V. CONCLUSIONS 

We have shown that theories of TeV-scale physics can 
have dark matter candidates whose annihilation proceeds 
via the formation of near-threshold WIMPonium bound 
states. Depending upon the closeness-to-threshold of 
the weakest-bound state, these can lead to a substan- 
tial velocity-dependent amplification of the dark matter 
annihilation cross-section, preferring the lowest velocity 
dispersion environments all other factors being equal, and 
providing a new dynamical source of "boost factors" for 
indirect detection signals. In addition, the amplified ra- 
diative capture to more deeply bound states, and the 
transitions among such bound states can both lead to a 
rich spectrum of discrete 7 lines which, if observed, would 
give striking confirmation of the mechanism. 

During the preparation of this work, [11] appeared. 
This paper also considers aspects of the phenomenology 
of WIMP bound states in the context of their annihila- 
tions and the consequences for dark matter indirect de- 
tection. 
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